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ABSTRACT 

Long-term JHK light curves have recently become available for large numbers of the 
more luminous stars in the SMC. We have used these JHK light curves, along with 
OGLE V and / light curves, to examine the variability of a sample of luminous red 
giants in the SMC which show prominent long secondary periods (LSPs). The origin of 
the LSPs is currently unknown. In oxygen-rich stars, we found that while most broad 
band colours (e.g. V—I) get redder when an oxygen-rich star dims during its LSP cycle, 
the J-K colour barely changes and sometimes becomes bluer. We interpret the J-K 
colour changes as being due to increasing water vapour absorption during declining 
light caused by the development a layer of dense cool gas above the photosphere. This 
result and previous observations which indicate the development of a chromosphere 
between minimum to maximum light suggest that the LSP phenomenon is associated 
with the ejection of matter from the stellar photosphere near the beginning of light 
decline. We explore the possibility that broadband light variations from the optical 
to the near-IR regions can be explained by either dust absorption by ejected matter 
or large spots on a rotating stellar surface. However, neither model is capable of 
explaining the observed light variations in a variety of colour-magnitude diagrams. We 
conclude that some other mechanism is responsible for the light variations associated 
with LSPs in red giants. 
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1 INTRODUCTION 


All luminous red giants appear to be variables and they 
are generally known as Long Period Variables(LPVs). They 
vary with periods which fall on at least eight period- 
luminosity sequenc es which have been labeled A', A, B , 
C', C, D, E and F (IWood et al.ll 19991:ISoszviiski et al.ll2004 
llta et al]l2004l : pTabur et al. 2010l -' S^zviiski fc Woodll2013ll . 
Sequences A'-C and F are thought to consist of stars pulsat¬ 
ing in the radial fundamental and overto ne modes and the 
non-radial dipo l e and quadrupole mode s (|Wood__eLaL 199^ 
Ita^et^ 2004 Talayam a et al.l l2013l : ISoszvhski k, WoodI 


20131: IStello et al.l |2014^ . These stars have been classified 


as Miras, Semi-Regulars (SRs) and OGLE S mall Amplitude 
Red Giants fOSARGsl l|Wood et "^119991 : ISoszvhski et al.l 
I2OO4I : ISoszvhski et al.ll2007l l. On the other hand, sequence E 
is known to consist of close binary system showing ellipsoidal 
and eclipsing variabi lity JWood et al.l 1 19991 : IPerekas et al.l 
I2OO6I : l^szvilskil [20071 1 . 

The last sequence, called sequence D, exhibits the 
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longest periods of all the sequences. Sequence D stars show 
photometric variations with periods of ~400-1500 days and 
they also show another variations with a shorter period that 
usually correspond to the period of sequence B. In these 
stars, the shorter period is called the primary period and 
the longer period corresponding to sequence D is known 
as a Long Secondary Period (LSP). LSP stars have been 
found within the LMG, SMC and Galactic populations of 
LPVs. In fact, approximately 25-50% of luminous red giants 
have been found to show LSPs in MACHO and OGLE dat a 
(|Wood et al.lll99!j : ISoszvhski et ahllioOTl : iFraser et al.ll200^ 1 . 

However, the origin of the LSP phenomenon is still unknown 
and many explanations have been suggested. Some of these 
are now highlighted. 

The most straightforward explanation is pulsation. 
However, all pulsation explanations have problems. Since, 
at the given luminosity, the typical LSP is approximately 
4 times longer than the period of a Mira (on sequence C), 
and since the Miras are known to be pulsating in the radial 
fundamental mode, an explanation in terms of normal-mode 
radial pulsation is not possible. Non-radial g-mode pulsation 
was proposed as the explanation of the LSP phenomenon 
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by IWood et alj (Il999l l since g-mode periods can be longer 
than the radial fundamental mode period. However, because 
red giants have a large convective envelope and only a thin 
outer radiative layer where g-modes could develop, the am¬ 
plitude of g-modes should be too small to correspond to the 
observed light and radial velocity amplitudes of LSP stars 
dWood et ahlliool i. 

Some researchers have also argued for a close binary hy- 
pothesis ( li ke seq uence E) to explain the LSP phenomenon. 
ISoszvrisM (l2007h found at least 5% of LSP stars showed 
ellipsoidal or eclipse -like variability. He also noted, as did 
iDerekas et al.l (l2006l l , that if the sequence E ellipsoidal vari¬ 
able were plotted using the orbital period rather than the 
photometric period (half the orbital period) then sequences 
D and E overlapped on the period-luminosity diagram. How¬ 
ever, the large amplitude of many LSP light curves (up to 
1 magnitude) would require a close pair of eclipsing red gi¬ 
ants, very unlikely in itself, and this would lead to alter¬ 
nating deep an d shallow minima w hich are not seen. On 
the other hand, IWood et al.l lll999ll found that many light 
curves due to the LSPs show faster decline and slower rise 
and they proposed that a binary with a companion sur¬ 
rounded by a comet-shaped cloud of gas and dust could give 
rise to those light curve features. Various observational and 
nu merical results have conflicted with the binary hypothe¬ 
sis. I Hinkle et al.l (l2002ll observed radial velocity variations 
in six local red giants with LSPs and found that the ve¬ 
locity curves all had a similar shape. For eccentric binaries 
with a random longitude of periastron, the light curve shape 
sho uld vary from star to star. The lack of such variability 
led iHink le et al. (l2002h t o rej ect the binary expla nation for 
LSPs. lWood et ahT i 2004h and iNicholls et al] ll2009lj obtained 
radial velocity curves for stars in the Galaxy and LMC and 
they argued that the small typical full velocity amplitude of 
~3.5 km s“^ meant that all binary companions would need 
a mass near 0.09 M©. Since companions of this mass are rare 
around main-sequence stars, it seems unlikely that the 25- 
50% of red giants which show LSPs could have a companion 
of this mass. 

Star spots due to magnetic activity o n AGB stars have 
been proposed by some researchers (e.g. ISoker fc ClavtonI 
Il999l j and a rotating spotted star has been suggested as 
an explanation of the LSPs. Possible evidence for mag¬ 
netic activity in stars with LSPs is provided by the ex¬ 
istence of Hg absorpt i on lin es whose strength varies with 
the LSP (|Wbod_et_aL l2004 l. However, from observations, 
lOlivier fc Wood ( 200f ) found that the rotation velocity of 
LSP stars is typically less than 3 km/s and they noted that 
this rotation velocity is too small to explain the periods of 
less than 1000 days which are often observed for the LSP. 
For example, a 3 km/s rotation velocity gives rise to the 
period of 2868 days for an AGB star with a ~170 Rq. A 
way to overcome this object ion would be to ha ve multiple 
star spots. In another study. IWood et ^ (l2004| j found that 
with a simple blackbody model, it was difficult for a rotat¬ 
ing spotted star to simultaneously reproduce both the am¬ 
plitude and colour variation observed in the MAGHO obser¬ 
vations of stars with LSPs. The use of model atmospheres 
rather than blackbodies, and allowing for limb darkening, 
could alter this conclusion. 

The giant convective cell model i s a relatively n ew ex¬ 
planation of the LSP phenomenon dStotherj l2010lj . This 


model is somewhat similar to the rotating spot model in 
that it involves variable amounts of the visible stellar sur¬ 
face having a temperature cooler than normal. The photo¬ 
metric variations seen in both the rotating spot model and 
the giant convective cell model should be similar. 

Another potential origin for the LSP is a variable 
amount of circumstellar dust. IWood et al.l (Il999l j proposed 
that the periodic dust formation found in th e circumstellar 
region s of t heoretical mo d els of AGB stars bv I Winters et al] 
(Il994h and iHofner et all lll995h could give rise to the light 
variation of the LSPs. Evid ence for dust aro u nd re d gi¬ 
ants with LSPs was found bv lWood fc Nicholld (l2009l j who 
showed that LSP stars tend to have a mid-infrared ex¬ 
cess. On near- to mid-infrared color-magnitude diagrams, 
the LSP stars lie in a region similar to that occupied by R 
Coronae Borealis stars which are known to be surrounded by 
patchy circumstellar dust clouds. This result suggests that 
the variation of LSP stars might be due to the ejection of 
dust shells that may be patchy. Similarly, in studying the 
effect of dus t eject ion on MACHO colours and magnitudes, 
IWood et al.l ll2004 l noted that patchy dust clouds seemed 
to be required in order to reproduce the observed colour- 
magnitude variations. 

In this paper, we utilize new long-term near-infrared 
JHKs light curves of stars in the SMC obtained with the 
SIRI US camera on the IRSF 1.4 m telescope lllta et al.l 
l2015lj . We combine the JHKs observations with V and 
I monitoring of SMC stars by the Optical Gravi tational 
Lensing Experiment (OGLE) (ISoszvriski et al.ll201l| j. Using 
this combined data, we examine light and colour variations 
of LSP stars over a wide wavelength range. We then test 
models involving dust ejection or variable spots to see if 
they can reproduce the observed light and colour variations 
over the wid e wavelength rang e. For our computations, the 
DUSTY code dlvezic et al.lll999h is used to model the effects 
of eje cted dust and the 2012 ver sion of the Wilson-Devinney 
code llWilson fc Devinnevlll97ll l is used to model the effects 
on the light and colour of a spotted star. 


2 THE OBSERVATIONAL DATA 
2.1 The near-IR data 

A long-term multi-band near-infrared photometric survey 
for variable stars in the Large and Small Magellanic Clouds 
has been carried out at t he South Africa n Astronomical Ob¬ 
servatory at Sutherland (llta iliDliQTi). The SIRIUS cam¬ 
era attached to the IRSF 1.4 m telescope was used for this 
survey and more than 10 years of observations in the near-IR 
bands J(1.25^tm), i7(1.63/rm) and As(2.14^m) band were 
obtained. In this work, we select the SMC stars from the 
SIRIUS database. Variable and non-variable stars in an area 
of 1 square degree in the central part of the SMC have been 
observed about 99-126 times in the period 2001-2012 and a 
total of 340147, 301841 and 215463 sources were identified 
in J, 77 and Ks, respectively. We note that the photometric 
detection range of the SIRIUS camera is about 8~18 mag¬ 
nitudes in Ks band. A total of 12008 variable sources of all 
kinds were detected of which 4533 w ere detected in a ll three 
wavebands (data release ver 130501. llta et al.l 1I2OI5IU . 
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2.2 The optical data 

We obtained V and I ba nd time series of SMC red giants 
from the OGLE project l|Soszvhski et al.ll201ll L Typically, 
about 1000 observing points were obtained in the I band by 
OGLE-II+III while a much smaller number of about 50-70 
poi nts were obtained in th e V band. 

ISoszvhski et ahl (l201lll divided the LPVs into LSP and 
non-LSP stars according to their positions in the period- 
luminosity relations for variable red giants. We found nearly 
700 LSP stars in the area of the IRSF infrared survey. In or¬ 
der to determine if an LSP star is oxygen-rich or carbon-rich, 
we initially adopted th e classification method introduced by 
ISoszvriski et al.[ (( 20 TJ). The stars we examined in detail had 
their oxygen-rich or carbon-rich status checked against op¬ 
tical spectral classifications in the literature. 

2.3 The combined near-IR and optical data 

By combining the SIRIUS data with the OGLE data, we 
obtained 298 oxygen-rich LSP stars and 88 carbon-rich LSP 
stars in the region monitored by the SIRIUS camera. For all 
sample stars, a small number of apparent bad data points 
were removed from each time series. Then, by using a first 
order Fourier fit to the I band time series, we obtained the 
period and amplitude corresponding to the largest ampli¬ 
tude mode of the light curves. In all cases, the obtained pe¬ 
riods and amplitu des were almost identi cal to those in the 
OGLE-III catalog (ISoszvhski et al.ll201llL Using the period 
from the fit to the I light curve, a fit was also made to the 
Ks light curve in order to derive its amplitude. 

Since our aim is to look for and to model colour and 
magnitude changes, we need the observed changes to be sig¬ 
nificantly larger than the noise in the data. The limiting 
factor in this study is noise in the near-IR data. We concen¬ 
trate on the brightest stars with Ks < 13 and for them the 
useful lower limit for clearly apparent Ks band variation 
is about 0.04 magnitudes. Similarly, we set a useful lower 
limit for I band variation of about 0.2 magnitudes. Figure[T] 
shows the full amplitude of the observed LSP stars in the / 
and Ks bands. This figure shows the selected stars, which 
having full amplitudes larger than 0.04 magnitudes in Ks 
and 0.2 magnitudes in I. This selection gave us a sample 
of 7 oxygen-rich stars and 14 carbon stars for analysis. The 
spectral types of our sample stars obtained from the OGLE 
photometry agree with the classifications in the SIMBAD as¬ 
tronomical database where the spectral types were obtained 
from spectroscopic determinations. The properties of these 
stars at maximum light are listed in Table [T] 

Near infrared JHKs photometric data as a function of 
Julian date (JD) for a sample of 7 oxygen-rich stars and 14 
carbon stars are published along with this paper. The file 
name of the time-series data indicates star-name and wave¬ 
band, like ”OGLE-SMC-LPV-07852.J.dat”. Table [5] shows 
time-series for a sample of the LSP star, which consists of 
seven columns as follows. 

Column 1 : Julius day. 

Column 2 : Calibrated (2MASS referenced) magnitude. 
Column 3 : Error in deferential image analysis in magnitude. 
Column 4 : Error in reference magnitude. 

Column 5 : Error in conversion offset in magnitude. 

Column 6 : Name of the survey region in which the variable 
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Figure 1. The relation between the I and Ks band full amplitude 
of 298 oxygen-rich stars and 88 carbon stars. The stars selected 
for detailed study have I amplitude greater than 0.2 magnitudes 
and Ks amplitude greater then 0.04 magnitudes. 


star is detected. 

Column 7 : Name of the survey subregion in which the vari¬ 
able star is detectecL _ _ 

for more detail, see llta et al. I (Eli). 


3 SEPARATING THE LSP VARIATION FROM 

THE PRIMARY PERIOD VARIATION 

Light curves of all LSP stars show variations with the pe¬ 
riod corresponding to sequence D and also variations with 
the primary period, which mainly corresponds to sequence 
B. Therefore, in order to study the LSP variation only, we 
need to separate each light curve into the components due 
to variation with the LSP and variation with the primary 
period. The separation method we used is now described. 

At first, we adopted a reference time to and then delim¬ 
ited the light curve into time intervals equal to the size of 
the LSP (P), with to as the start of an interval. Then we fit¬ 
ted the observed data points in each interval independently 
using a second order Fourier series Po with the period equal 
to that of the LSP, 

2 

Fo{t) = Ao 4- ^ ao sin ^27ri ^ p^° ) ’ 

i=l 

The fit was performed on independent intervals of length 
P rather than the full light curve interval because the light 
curves of LSP stars vary substantially from cycle to cycle of 
the LSP. Furthermore, in order to minimize the dependence 
of the final fit on the adopted value of to, the fit was made us¬ 
ing 4 values for to, tq, ri, T 2 and ra, with each of these values 
increasing by \P, giving fits Foit), Pi(t), F 2 {t) and F^it), 
respectively. In this study, tq corresponds to JD2450000. A 
requirement for a fit to be made was that there be at least 6 
points in each of the first and second halves of the fit inter¬ 
val. If such points were not available, no fit was made. The 
final adopted fit F(t) at the observed times is given by 
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Table 1. Basic parameters of the program stars at maximum light 


Name 

RA 

Dec 

V 

(mag) 

I 

(mag) 

J 

(mag) 

H 

(mag) 

Ks 

(mag) 

L 

(Lo) 

Tbb 

(K) 

Teff 

(K) 

Plsp 

(days) 

Sp. 

type 

OGLE-SMC-LPV-07852 

12.402540 

-72.97075 

16.60 

14.34 

12.91 

12.05 

11.73 

4617 

2803 

3434 

645 

O 

OGLE-SMC-LPV-09856 

13.186620 

-72.95178 

16.49 

14.37 

12.96 

12.10 

11.86 

4408 

2928 

3558 

615 

O 

OGLE-SMC-LPV-10774 

13.530120 

-72.47989 

- 

14.98 

13.73 

12.81 

12.57 

2170 

2844 

3474 

396 

o 

OGLE-SMC-LPV-12332 

14.166840 

-73.07389 

- 

14.26 

12.95 

12.08 

11.82 

4455 

2887 

3517 

632 

o 

OGLE-SMC-LPV-14045 

14.987085 

-72.94936 

17.50 

15.13 

13.72 

12.79 

12.57 

2197 

2865 

3495 

390 

o 

OGLE-SMC-LPV-14715 

15.332205 

-72.86369 

- 

14.27 

12.92 

12.03 

11.75 

4544 

2815 

3445 

601 

o 

OGLE-SMC-LPV-14781 

15.361920 

-72.85678 

16.48 

14.62 

13.30 

12.44 

12.23 

3245 

2991 

3620 

469 

o 

OGLE-SMC-LPV-08199 

12.555540 

-73.18697 

17.75 

14.90 

13.29 

12.27 

11.64 

3745 

2203 

2792 

767 

c 

OGLE-SMC-LPV-08280 

12.583995 

-72.80533 

16.36 

13.78 

12.33 

11.36 

10.93 

7922 

2490 

3106 

1008 

c 

OGLE-SMC-LPV-09341 

12.990495 

-73.20256 

- 

13.92 

12.48 

11.63 

11.26 

6419 

2734 

3362 

1081 

c 

OGLE-SMC-LPV-10109 

13.281915 

-73.12992 

16.08 

13.84 

12.52 

11.60 

11.20 

6411 

2592 

3215 

964 

c 

OGLE-SMC-LPV-12653 

14.317380 

-72.80128 

- 

14.83 

13.10 

12.03 

11.44 

4465 

2200 

2789 

1356 

c 

OGLE-SMC-LPV-13340 

14.654580 

-72.44625 

16.28 

14.12 

12.59 

11.66 

11.28 

6003 

2606 

3230 

1014 

c 

OGLE-SMC-LPV-13557 

14.755710 

-72.46908 

16.31 

13.91 

12.45 

11.54 

11.13 

6876 

2592 

3215 

1032 

c 

OGLE-SMC-LPV-13748 

14.847210 

-73.04900 

- 

14.69 

13.09 

12.11 

11.67 

3977 

2468 

3082 

925 

c 

OGLE-SMC-LPV-13802 

14.871165 

-72.65750 

17.08 

14.19 

12.69 

11.71 

11.16 

6085 

2334 

2936 

1011 

c 

OGLE-SMC-LPV-13945 

14.946045 

-72.71094 

- 

14.06 

12.63 

11.73 

11.29 

5872 

2560 

3180 

927 

c 

OGLE-SMC-LPV-14084 

15.004245 

-72.89839 

16.30 

14.27 

12.74 

11.83 

11.53 

5042 

2759 

3388 

799 

c 

OGLE-SMC-LPV-14159 

15.048705 

-72.66822 

- 

14.42 

12.96 

12.04 

11.57 

4406 

2495 

3112 

913 

c 

OGLE-SMC-LPV-14829 

15.382875 

-73.09919 

15.69 

13.70 

12.37 

11.54 

11.21 

6950 

2839 

3469 

780 

c 

OGLE-SMC-LPV-14912 

15,434040 

-73.29506 

- 

14.49 

12.92 

11.91 

11.31 

5135 

2244 

2837 

722 

c 


Notes: Tr b is the blackbody temperatur e derived from J — Kg while Tgff is the effective temperature derived from J — Kg using the 
formula in lBesselh Wood. &: Evan j lll983l) . Luminosities L were derived by applying a bolometric correction to Kg (see Section [4Jl . 


Table 2. The first ten records from the time series data for a variable source in J band (OGLE-SMC-LPV-07852). The seven columns 
are explained in the text. The full version of this table is available as online material. 


Column 1 
(day) 

Column 2 
(mag) 

Column 3 
(mag) 

Column 4 
(mag) 

Column 5 
(mag) 

Column 6 

Column 7 

2452092.635160 

12.930 

0.003 

0.019 

0.003 

SMG0050-7250 

G 

2452212.344921 

12.959 

0.002 

0.019 

0.003 

SMG0050-7250 

G 

2452213.302422 

12.958 

0.003 

0.019 

0.003 

SMG0050-7250 

G 

2452214.397972 

12.940 

0.002 

0.019 

0.003 

SMG0050-7250 

G 

2452246.289626 

12.991 

0.003 

0.019 

0.003 

SMG0050-7250 

G 

2452285.324709 

13.007 

0.003 

0.019 

0.003 

SMG0050-7250 

G 

2452298.305372 

12.974 

0.003 

0.019 

0.003 

SMG0050-7250 

G 

2452390.680728 

12.931 

0.002 

0.019 

0.003 

SMG0050-7250 

G 

2452411.689245 

12.948 

0.002 

0.019 

0.003 

SMG0050-7250 

G 

2452415.634797 

12.965 

0.003 

0.019 

0.003 

SMG0050-7250 

G 


3 

j=o 

When there were fewer than 4 fits, the average was taken 
over the nnmber of fits available. An example of the final 
light curve fits for each band in a representative star is shown 
in Figured 

As one can see from Figure[2l observation times of the 
OGLE I band usually differed from those of the OGLE V 
band and the SIRIUS J, H and Ks bands, the latter being 
taken simultaneously. The I band observations are always 
the most frequent. In our analysis, we will examine colours 
relative to the I band magnitude. In order to derive the 
colours and I band magnitudes at simultaneous dates, the / 
band fits to the LSP variations were used. Magnitudes and 
colours were derived from the fits at observation times of 


the V band and the Ks band. As noted above, fits are not 
available at all times in a given band. 


4 MODELING THE LSP MAGNITUDE AND 
COLOUR VARIATIONS 

We will now create models for the magnitude variations in 
various bands associated with LSPs in order to see if the 
assumed models are consistent with the observations. The 
two models we examine are obscuration by dust and the 
variation of spots on a star. 

4.1 The dust model 

One possible explanation of the LSP phenomenon is episodic 
circumstellar dust absorption. Here we assume that the cen- 
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HJ 0-2450000 

Figure 2. The complete light curves of OGLE-SMC-LPV-10774 
(P=396d) in the bands VIJHK. The raw data points are shown 
as small filled grey triangles while the fit values are shown as large 
filled black circles. 

tral star ejects mass for several hundred days in each LSP 
cycle, leading to the formation of thin spherical mass shells. 
Dust grains which form in the ejected matter will be carried 
along with the expanding mass shells leading to dimming 
of the light from the central star. We assume that the dust 
shells are ejected spherically although this is not an essential 
aspect of the dust model for the optical and near-infrared 
bands which we are examining since thermal emission in 
these bands is not strong. 

The amount of dimming of a star due to circumstellar 
dust depends on the optical depth through the dust. For a 
spherical dust shell with an outer radius i?out (corresponding 
to the start of a mass ejection episode), an inner radius i?in 
(corresponding to the end of a mass ejection episode, or 
the dust formation radius in a newly-forming inner shell), 
a constant expansion velocity v and a mass ejection rate 
M, the optical depth t\ throu gh the shell is given by (e.g. 
lOroenewegen fc de .londllQD.^ l 

M'i'Qx/a , 1 _ 1 ^ 

upg ^Rin Rout^’ 

where T is the dust to gas mass ratio, Qx is the absorption 
efficiency of the grain at wavelength A, a is the grain radius 
and pg is the grain density. _ 

We have used the dusty code dlvezic et al.l Il999l j to 
examine the effect of a spherical dust shell on the broad 
band colours of LSP stars, dusty computes the SED of a 
dust enshrouded LSP star for one given configuration of the 
dust shell. From Equation [T] we can see that the innermost 
dust shell plays the most important role for dust extinction 
when multiple dust shells of similar origin exist. In our model 
we only consider the innermost dust shell. 

The first input parameter for dusty which we consider 
is the dust shell thickness. AGB stars are well-known mass¬ 
losing stars in whic h the terminal vel o city of the mass flow 
is ~ 10-20 km s ~^ JPecin et al. I l2007l : iDe Beck et Zll20ld: 
iLombaert et al.ll2013l j. On the other hand, R Coronae Bore¬ 
alis (RGB) stars are similarly luminous evolved sta rs which 
eject mass semi-periodically in random directions dClavtorJ 
Il99d j. The large-amplitude variability of these stars is due 


to dust absorption by those ejected dust clouds which lie in 
the line of sight to the star. Observation ally, the ejection ve¬ 
locity of the clouds is '^200-250 km s“^ (lFeasfj|l97^ : IClavtorJ 
ll996l:lFeastlll990lL In view of these velocity estimates, and 
since we do not know what causes dust shell ejections in LSP 
stars, we considered two values for the wind velocity v of 15 
and 220 km s“^, and the velocity is assumed to be constant 
with radius. We assume that dust is ejected for a time inter¬ 
val equal to half the length of the LSP so the shell thickness 
is uP/2, where P is the length of the LSP. Our models are 
thus made with two values for the shell thickness. 

The second input parameter we require for dusty is 
the temperature Tin at the inner edge of the dust shell. 
We assume Tin is the dust condensation temperature. The 
dust condensation temperature is c ommonly assumed to 
lie in the ran ge of 800-1500 K (e.g. Ivan Loon et al.l l2005l : 
ICassara et al.ll2013l l and in our models we consider Tin val¬ 
ues of 800 and 1500 K. We note that Tin determines the 
inner radius of the dust shell. 

The third input parameter for dusty is the grain type. 
The assumed chemical composition of the grains is dif- 
ferent for mod e ls of oxygen-rich stars and carbon stars. 
Ivan Loon et al.l (l2005l ) used dusty to empirically determine 
possible grain types for oxygen-rich and carbon AGB stars 
in the LMG by requiring an acceptable fit to the SED of the 
observed stars. We use the grain chemical composition that 
they adopted. For carbo n stars, we us e dust grains of 50% 
amorphous ca r bon f rom iHanne 3 (Il988lj . 40 % graphi t e from 
iDrairie fc Led (Il984l l and 10 % SiC from IPegourid (Il988l j 
a nd for oxygen-r i ch sta rs, we use the astronomical silicates 
of iDraine fc Led (1 19841 1. 

The stellar luminosity L and the effective tempera¬ 
ture Tefi of our LSP stars are also required as input to 
DUSTY. These were computed from the J and Ks magni¬ 
tudes at maximum light, when circumstellar extinction is a 
minimum. The Tefr values for the observed stars were es¬ 
timated from the J-Ks color at maximum light using the 
formula given in iBessell. Wood, fc Evans! (Il983l l. We also 
used J-Ks to calculate a blackbody temperature Tbb for 
the central star at maximum. Luminosities were computed 
using a bolometric correction BCk to the Ks magnitude. 
For oxygen-rich stars, we adopted the formula BCk ~ 
0.60-1-2.56(J—A)—0.67(J—A)^ given bv iKerschbaum et al.l 
1 I 2 OIOII and for carbon stars we adopt the formula BCk = 
1.70 -f 1.35( J — K) — 0.30( J — A)^ given bv IBessell fc Woo'd! 
(Il984ll . To estimate of absolute bolometric magnitude, we 
used the distance mo dulus to the SMC of 18.93 given by 
iKeller fc WoodI (l2006fl . The derived values for L, Tbb and 
Tefi are listed in Table [T] When using dusty, we assumed 
that the central star emits a blackbody spectrum of temper¬ 
ature Tbb . We are only interested in changes in magnitudes 
and colours, rather than the absolute values of the colours 
and magnitudes, and the blackbody assumption will have 
only an insignificant effect on these changes. 

The final input parameter to dusty is the optical depth 
Tv in the V band. The output from dusty is the SED of 
the central star extincted by a shell of visible optical depth 
Tv- Unless there is significant emission from the dust shell 
in each band under consideration, for a given value of tv 
the extinction of the central star in each band is essentially 
independent of the details of the shell radius or thickness. 
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Figure 3. Comparisons of observations with models of expansion velocity 15 km s~^. The 7 panels of the left column represent the I 
and I — J variations for the oxygen-rich stars and the 14 panels of the middle and right columns are for the carbon stars. The red points 
show the LSP variations obtained from Fourier fits to the observed light curves. The blue lines show least squares fits to these points 
(the line is omitted if the relative error in the slope is more than 0.2). The black continuous lines with open circles show the magnitude 
and colour variations for models with an inner dust shell temperature of 800 K while the triangles are for models with an inner dust shell 
temperature of 1500 K. The top end of each curve corresponds to shell optical depth rv = 0 and the points are spaced at intervals of 
0.5 in TV = 0. The black dotted lines show the locus of un-extincted blackbodies of different temperatures and the constant luminosity 
given in Tabled for each star. Note that the scale of the vertical axis is the same for all oxygen-rich stars and for all carbon stars. 
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Figure 4. The same as Figure [3| but for models with an expansion velocity of 220 km 
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4.2 Comparison of the dust model with 
observations 

4-2.1 The variation of I with I — J 

In general, the I light curves are the best sampled of the 
light curves in the different filters, the J, H and K light 
curves are the next best sampled and the V light curves 
are mostly only sparsely sampled. Thus when comparing 
observations against models in magnitude-colour plots, we 
use I as the magnitude and a colour involving I and another 
magnitude. The other magnitude should be one of J, H or 
K because of the higher frequency of observation compared 
to V. Ideally, I — K would be best because of the large 
wavelength range spanned by this colour, but as we will 
show below, the K band in O-rich stars is contaminated by 
variable H 2 O absorption, as is the H band, so we use I — J 
as the colour in our first comparison of observations and 
models. 

In Figures|3] and |4] we show the I magnitude plotted 
against the I — J color for the observed samples of oxygen- 
rich and carbon stars. Figure^ displays models with a shell 
expansion velocity of 15 km s“^ (thin shell) while Figure|4] 
displays models with a shell expansion velocity of 220 km 
s“^ (thick shell). Note that it is not the absolute values of 
points in the diagrams that is important, it is the slope of 
the line representing the variation of I with I — J that mat¬ 
ters. If the LSP variations are caused by variable amounts of 
circumstellar dust alone, then the slope will depend almost 
entirely on the dust properties and not on the spectrum of 
the central star. The central star, which in these models we 
assume to emit as a blackbody with the L and Tbb given in 
Table [1] determines the position of the model in Figures[3] 
and m when there is no circumstellar extinction. 

Firstly, we make some comments on the models. In both 
Figure|3]and Figure|4j the two model lines differ slightly due 
to re-emission from the dust shell at J and to a lesser extent 
I. Comparing Figure|3]with Figurejd] corresponding to mod¬ 
els with different dust shell thickness, we see that there is a 
slightly faster increase in I — J with I for the thinner shell 
of FigureO especially for the hotter inner dust shell temper¬ 
ature of 1500 K. The reason for this is that for the thinner 
dust shell there is more emission at J. Also, for a given op¬ 
tical depth Tv to the centre of the disk of the observed star, 
the total extinction averaged over the stellar disk is slightly 
greater for the thinner shell in Figure[3]than for the thicker 
shell in Figure|4](the lines in Figure[3]are longer than in Fig¬ 
ure 0. This is because the extinction near the limb of the 
disk is smaller relative to the extinction at the centre of the 
disk for a thicker shell. 

The observed variation of I with I — J associated with 
the LSP (red points) is generally consistent with the slope of 
one of the model lines in Figures|3]or0 This indicates that 
the variation of I and J due to an LSP could be explained 
by variable dust absorption alone. 


4-2.2 The variation of I with J — Ks 

Figures[5] and [6] show the variation of I with J — Ks- As 
one can easily see, the observational variations due to the 
LSPs are inconsistent with the model variations in all cases. 
Most importantly, for the O-rich stars, the observed J — Ks 
colour tends to get bluer as the star gets fainter in I. There 


is no way that such variations can be the result of variable 
amounts of dust extinction only. 

The most probable explanation for the variation of 
J — Ks in O-rich LSP stars is strong H 2 O absorption in 
Kg band. Significant broad H 2 O band absorption is found 
between about 1.4^m and on the edges of the H 

and Ks filter bandpasses. In non-variable stars, this ab¬ 
sorption is seen f rom about M6 HI to later spec tral type 
in M giants (e.g. iTsuiil [ 2 O 0 O : iRavner et al.l [200^ . Strong 
water absorption in the H and K bands is also com¬ 
monly seen in l arge amplitude pulsating red giants (Mira 
variables) (e.g. I Johnson fc Mended 1 1970l i. The reason for 
this absorption is the existence of a relatively dense shell 
with a temperature of about 1800 K elevated above the 

S ihotosphere by shock waves as sociated with the pulsation 
Ireland. Scholz. fc_Wood 2011 ). Using the models atmo¬ 
sphere results of Houdashelt et all (l2000li . we find that the 
V — I colours of the stars we are examining here yield 
Teff ~3400-3600 K, corresponding to spectral types M3“M5. 
Due to these relatively early M spectral types, the results 
noted above mean that strong water absorption in the Ks 
band can not occur in a hydrostatic atmosphere. The ele¬ 
vation of a shell of mass above the photosphere is required 
to produce strong water absorption, as in the case of Mira 
variables. We suggest that this shell elevation, which may 
also lead to dust formation, occurs at the beginning of the 
luminosity declines associated with the LSPs, leading to the 
formation of increasing H 2 O absorption in the Kg band and 
blueing of the J — Kg colour. 

This result is a second piece of observational evidence 
for the existence of a disturban c e abo ve the photosphere 
caused by the LSP. IWood et al.l (l2004h found variable Ha 
absorption in stars with LSPs which they attribute to a chro¬ 
mosphere of variable strength, with the strongest Ha absorp¬ 
tion near maximum light. Combined with the results above 
for the H 2 O absorption layer, we have a picture whereby 
matter is ejected above the photosphere beginning near the 
luminosity decline of the LSP forming a relatively dense shell 
of temperature ~1800K. As the LSP cycle progresses, this 
shell appears to be slowly heated by non-thermal processes 
to ~8000 K giving rise to a chromosphere which is then dis¬ 
placed by the next ejected shell containing H 2 O. 

Turning to the carbon stars, we see that they do get 
redder as they get fainter in I but in all cases they clearly do 
not redden as much as predicted by the models. To bring the 
observed variation in J —Kg with decline in I into agreement 
with the models, while keeping the variation of / — J with I 
in agreement with the models, would require a strange dust 
opacity with a large bump in the Kg band. We do not know 
of any dust composition that could produce such a bump 
(see Section 14.2.41 for a discussion of the effect of different 
grain opacities). 

We note that ejected clouds covering only part of the 
stellar surface rather than the complete surface as for spher¬ 
ical shells can not explain the observed variation in J — Ks 
with declining I. For all but very large optical depths, the 
J — Kg colour in a cloud model increases slightly faster with 
decline in I than for a spherical shell model, whereas the ob¬ 
servations require a slower increase in J — Kg colour. Also, 
given that the variation of the observed I — J colour with 
decline in I agrees with the spherical shell models, a cloud 
model will not agree with the observations (the V — I colour 
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Figure 5. The same as Figure[3]but using the J — Ks colour rather than I — J. 
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Figure 6. The same as Figure[4|but using the J — Ks colour rather than I — J. 
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in the cloud model increases more slowly with decline in I 
than in the spherical shell model). 

In summary, our conclusion from examining the varia¬ 
tions of LSP stars in the {I,J — K) plane is that, for both 
0-rich and C-rich stars, a model involving only variable dust 
extinction is not viable as an explanation for LSP behaviour. 


4 ..2.3 The variation of I with V — I 

Although the V observations are not as frequent as observa¬ 
tions in other bands, we show the variation of / with V — I 
for completeness in Figures [7] and [S] For both the V and I 
bands, re-emission of light absorbed by the shell is weak so 
the model lines corresponding to the two inner shell tem¬ 
peratures of 800 K and 1500 K are close together in both 
figures. 

For most 0-rich stars, the observed V — I tend to in¬ 
crease slightly more rapidly with declining I than in the 
models. For most C-rich stars, the observed and model varia¬ 
tions have consistent slopes although for OGLE-SMC-LPV- 
13557 and OGLE-SMC-LPV-13945 the slopes are different. 
Overall, the agreement between the observations and dust 
model predictions is not good in the {I,V — I) plane but the 
results are not dehnitive. 


4-2.4 Testing different dust grains 

Different dust grains have different absorption and scatter¬ 
ing coefficients. Thus we tried changing the dust grain com¬ 
position to see if this would give agreement between the 
light variations in the dust shell model and the observed 
light variations, especially in the I and J — Ks plane for the 
carbon stars. We examined four dust grain types available 
in DUSTY (see Table O for oxygen-rich stars and for carbon 
stars, including the grain types we used in Section [4.II The 
thin dust shell model with an expansion velocity of 15 km 
s“^ was adopted. 

Figure [9] shows comparisons of observed variations in 
the I and J — it's diagram with dust shell models for the 
different grain types. The model slopes for oxygen-rich stars 
can be changed appreciably by different dust grains but none 
of the models reproduce the observed slopes. In particular, 
none of the models can explain the observation that the 
J — Ks colour get bluer as the star gets fainter in I. Turning 
to the carbon stars in Figure [51 we see that changing the 
grain type in the models changes the slope of the variations 
but again none of the adopted grain types can consistently 
reproduce the observed slope. These results strengthen the 
hnding of Section 14.2.21 that dust shell models can not by 
themselves explain the broadband light and colour variations 
in LSP stars. 


4.3 The dark spot model 

Another possible explanation of the LSP phenomenon is 
a dark spot on a rotating AGE star. ISoszvhski fc Udalskil 
1 I 2 OI 4 II noted that the observed light curves of LSP stars are 
similar to those of the known spotted stars. Moreover, the V 
to I amplitude ratios of LSP stars and known spotted stars 
are similar and they are different from the ratios seen in 
the ellipsoidal binaries of sequence E which lie near the LSP 


Table 3. Dust composition 


Sp. type model name 

composition 

0 rich star 

model 0 

100% ^astronomical silicate 


model 1 100% ^warm 0-deficient silicate 


model 2 

100% ^FeO 


model 3 

100% '^glassy olivine 

C rich star 

model 0 

50% ^amorphous carbon, 

40% ^graphite, 10% ®a-SiC 


model 1 

100% ^graphite 


model 2 

100% ^amorphous carbon 


model 3 

100% ®o-SiC 

Note: toraine & Led lll984h. 

dOssenkoof et al.l lll992ri. 

^Hennine: et alj i 

1995^. toorschner et alj lll995^. ^Hanned 


stars in period-luminosity plots. Here we test a simple rotat¬ 
ing dark spot model for LS P stars using the latest ver sion of 
the Wilson-Devinney code l| Wilson &: Devinnevll97l]l which 
has the ability to simulate the light curve of a spotted star 
in a binary system. Since our aim is to examine the light 
and colour variations due to the appearance of a dark spot 
on a red giant, our binary systems consist of a rotating red 
giant star with a circular isothermal spot and a very faint 
orbiting companion which can not be seen. Hence, we obtain 
the light and colour variations due to only the dark spot on 
the rotating red giant. For the red giant star, we assume 
a mass of ~1.0Mq and a radius of ~17O7?0, which gives a 
surface gravity \ogg [cgs] of —0.02. A metallicity, log[ Fe/H] 
= —1.0 and a gravity-darkening coefficient, gi=0.32 (ILucvI 
are used. For the calculation o f oxygen-rich stars, we 
use th e model stellar atmospheres of Ivan Hamme fc Wilson! 
1 I 2 OO 3 II while for carbon stars, a blackbody is used. We note 
that the model atmospheres do not take into account the 
effect of the absorption due to H 2 O molecules. 

The amplitudes of the light variations due to the pres¬ 
ence of spots are determined by the spot sizes and the 
temperature differences between the spots and the normal 
stellar surface. The temperatures of spots on AGE stars 
are not kn own well. For the spot temperature, we follow 
the work of ISoker fc GlavtonI (Il999ll . They initially adopted 
the solar value of | for the ratio of spot to normal photo- 
spheric temperature. For an assumed photospheric tempera- 
ture of 3000 K, the spot temperature would then be 2000 K. 
ISoker fc Glavto^ (Il999l l then considered a range of spot tem¬ 
peratures from 1600 K (cool spot,) and 2500 K (warm spot) 
corresponding to a ratio of spot to normal photospheric tem¬ 
perature of 0.533 and 0.833, respectively. To cover a range of 
spot temperatures, we also make models with these two tem¬ 
perature ratios. We use these models to simulate each of the 
stars in Table [T] and we use the effective temperatures T^g 
in Table [T] as the temperature of the normal photospheric 
regions. 

The spot sizes for the two temperature cases are fixed by 
requiring that the amplitude in the I band is similar to the 
observed largest amplitudes, which are ~0.5 magnitudes for 
the oxygen-rich stars and >^0.7 magnitudes for the carbon 
stars. Table [4] gives the parameters of our spot models. It 
has been assumed that the inclination of the model rotation 
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Figure 7. The same as Figure [3] but using the V — I colour rather than I — J. 
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Figure 8. The same as Figure [4| but using the V — I colour rather than I — J. 
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Figure 9. The same as Figure[5]but including different dust compositions. For both oxygen-rich and carbon stars, the black solid lines 
with open circles and filled triangles correspond to the original dust compositions (model 0 in Table [Sj), while the green, magenta and 
blue colours correspond to models 1, 2 and 3, respectively. 
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Table 4. Spot temperatures and the spot sizes 


Sp. type 

temperature ratio 

spot radius(°) 

model 

0 rich star 

0.533 

35 

cool 


0.833 

45 

warm 


0.913 

90 

hemisphere 

C rich star 

0.533 

40 

cool 


0.833 

50 

warm 


0.888 

90 

hemisphere 


Note: Spot radius is defined as the angle corresponding to the spot 
radius when viewed from the stellar center. 


axis is 90 degrees to the line of sight and the center of the 
spot is on the equator of the star. 

The spot model also effectively simulates light varia¬ 
tions caus ed by the gian t convective cells which were pro¬ 
posed by llStothersll2010li to explain LSPs in AGB stars. In 
this model, the light variations are caused by the turnover 
of cooler and hotter giant convective cells at the stellar sur¬ 
face. To simulate this situation, we examine a model where 
the rotating dark spot covers a full hemisphere of the star 
(hereinafter, the hemisphere model). To reproduce the ob¬ 
served largest amplitudes in I band, we use the temperature 
ratios of the spot area to the normal area given by Table ID 
As in the case of the other spot models, we use the effec¬ 
tive temperatures Tefj given in Table [T] for the temperature 
of the normal photospheric regions in both oxygen-rich and 
carbon stars. 

4.4 Comparison of the spot model with 
observations 

4-4-i The variation of I with colours 

In Figures fTOl 1111 and 1121 we show the I magnitude plotted 
against the colours I — J, J — Ka and V — I, respectively, 
for the observed samples of oxygen-rich and carbon stars. 
The magnitudes of models at the maximum light, which 
corresponds to when there is no visible spot, are adjusted in 
the plots to the magnitudes at the observed maximum light 
for each star. 

As with our comparison of dust shell models with obser¬ 
vations, we test the spot models by looking for differences in 
the slopes in (/,colour) diagrams of model lines and obser¬ 
vational variations. Models with different ratios of the spot 
to normal photospheric temperature show different slopes. 
The models with higher spot temperature (i.e. larger tem¬ 
perature ratio) show shallower slopes. 

For carbon stars, the slope of the line representing the 
observations generally shows a good agreement with the 
model lines in each of Figures (TO] [11] and [12] (two exceptions 
are OGLE-SMC-LPV-13557 and OGLE-SMG-LPV-13945 in 
Figure [12] which show disagreement with their models). Ac¬ 
cording to the figures, the models with hotter spots (tem¬ 
perature ratio 0.833) and the hemisphere models (which also 
have a temperature ratio of 0.888) fit the models better in 
each of the (7,7 — J), (7,J — Ks) and {I,V — 7) diagrams. If 
sequence D light variations are due to rotating spots, these 
results suggest that the spot temperatures in carbon stars 
should be roughly 2500 K. 



I-J 

Figure 13. The variation of 7 with 7 — J in hemisphere spot 
models for oxygen rich stars when different inclination angles tor 
the rotation axis are assumed. 


In distinct contrast to the carbon stars, the oxgen rich 
star variations do not agree with model variations in any of 
the (7,7 — J), (7,J — Kb) or {I,V — I) diagrams. Although 
agreement would not be expected in the (7,J — Ka) plane 
where the observations are affected by absorptions due to 
H 2 O molecules, the bad disagreement in the (7,7 — J) and 
(7,1/—7) diagrams indicates that spot models cannot explain 
the light and colour variations due to the LSP phenomenon 
in oxygen-rich stars. Since the LSPs of both oxygen-rich 
stars and carbon stars should have the same origin, the 
failure of the spot models to reproduce the observations in 
oxygen-rich stars means that spots are not a reasonable ex¬ 
planation for LSPs in sequence D stars . This result agrees 
with the conclusion of lOlivier fc Wo^ (l2003h that the ro¬ 
tation periods of LSP stars are too long for a spotted star 
model to explain LSP variations. 


4.4-2 The effect on spot models of the inclination of the 
rotation axis and the spot lattitude 

The models above were made assuming spots on the equa¬ 
tor of a rotating star whose rotation axis is perpendicular 
to the line of sight. Changing the lattitude of a spot will re¬ 
duce its apparent area and hence the amplitude of the light 
variation it will cause. However, the slope of the locus of 
colour-magnitude variation should be unchanged since the 
reduction in apparent area also occurs as a spot on the equa¬ 
tor rotates towards the limb of the star (the very minor ef¬ 
fect of limb darkening will produce very small differences 
between the two cases). The same argument applies in the 
case of a star with a tilted rotation axis. We thus expect the 
results of the previous section to apply regardless of the in¬ 
clination of the rotation axis and the lattitude of the spots. 
As a test calculation, we examined the change in the 7 and 
I — J magnitudes for a star with the rotation axis tilted at 
30, 60 and 90 degrees. The results are shown in Figure [Till 
It can be seen the the slope of the variation is the same in 
all cases. 
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Figure 10. Comparisons of observations with dark spot models in the (7,7 — J) plane. The 7 panels of the left column show oxygen-rich 
stars and the 14 panels of middle and right columns show carbon stars. Red points and blue lines represent observations, as in Figure^ 
The black lines with open circles and the lines with filled triangles show the magnitude and colour variations for models with cool and 
warm spot temperatures, respectively, while the lines with crosses are for models with a hemisphere spot. 
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Figure 11. The same as Figure FTOl but using the J — Kg colour rather than I — J. 
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Figure 12. The same as Figure fTOl but using the V — 1 colour rather than I — J. 
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4.5 A combined dust and dark spot model 

In Section the possibility that a dust shell model alone 
could explain the LSP phenomenon was examined and re¬ 
jected. Similarly, in Section|4]4]a rotating spotted star model 
was examined and rejected. Here we examine the possibil¬ 
ity that a combination of a rotating spotted star and dust 
absorption at certain phases of the rotation could cause the 
observed behaviour of stars with LSPs. 

In this combined model, the effective temperature of 
the central star is assumed to vary sinusoidally to simulate 
a rotating spotted star and the optical depth of the dust 
shell is assumed to vary sinusoidally to simulate periodic 
dust absorption. The following equations are used: 

Tes{(l>) = Tefr,o + ATeff sin (27r(^-I-/3) (2) 

Tv {4>) = Tv,o sin {ncj)) (3) 

where ATefi and tv,o are the amplitudes of variation of the 
effective temperature and the optical depth in V of the dust 
shell, respectively. Teff.o is the average effective temperature. 
As a direct consequence of these equations, the central star 
luminosity is given by 

m = (4) 

where Lq is the average luminosity. The fluxes emitted by 
the central star in the various photometric bands were ob¬ 
tained assuming the central star was a blackbody. The dust 
absorption in photometric bands other than V was obtained 
using DUSTY. The dust compositions used in Section 14. II for 
oxygen-rich and carbon-rich stars was adopted and the shell 
outer radius was set to 1.5 times the inner radius, which 
roughly corresponds to the models described above with an 
expansion velocity of 15 km s“^. The inner dust shell tem¬ 
perature was set to 800 K. 

In the combined model, 4> varies from 0 to 1 as the 
spotted star makes one complete rotation. The phase shift /3 
between the rotation phase and dust absorption phase has 
been introduced to allow for dust absorption at different 
possible rotation phases. Since we do not know how dust 
absorption might be linked to the rotation phase, models 
were made with different values for j3. Note that the dust 
absorption goes through one maximum as the star makes a 
full rotation. 

Table [5] gives the parameters considered for the models 
described here. Two values of ATefi, 50K and 100K, were 
used for both oxygen-rich and carbon-rich stars. The values 
of Tv,o were adjusted so that the variations in magnitude 
caused by dust dominated for the models with ATefi = 50 K 
while the variations in magnitude caused by Teg variations 
dominated when ATefi = 100 K. 

Figure [2] shows the I magnitude plotted against V — I, 
I — J and J — K for the various models computed here. In 
general, the slopes of the variations in the (7,colour) dia¬ 
grams for the combined models lie between the slopes for 
pure dust absorption and pure Teg variations (spot model). 
The combined models do not generally represent an improve¬ 
ment over the pure dust or spot models. In particular, they 
can not explain the observation that J — K in oxygen-rich 
stars get bluer as the star dims. 


Table 5. and tv,o for models 


Sp. type 

7;fi.o(A) 

ATefi (K) 

TV,0 

model 

0 rich star 

2800 

100 

1 

modelO 


2800 

50 

2 

modell 

C rich star 

2500 

100 

0.5 

model2 


2500 

50 

1 

models 


5 SUMMARY 

We have examined broadband observations of samples of 
oxygen-rich and carbon stars in the SMC which exhibit LSP 
light variations. The observations are in the VIJHKs bands 
and they were obtained by the OGLE project and by the 
SIRIUS camera of the IRSF 1.4m telescope. 

The J and Ks observations of oxygen-rich stars reveal 
a significant new feature of LSP light variations. As light 
declines, the J — Ka colour barely changes or sometimes 
gets bluer. We interpret this as indicating that cool gas con¬ 
taining a significant amount of H 2 O is levitated above the 
photosphere of LSP stars at the beginning of light decline. 
Absorption of light in the Ka band by H 2 O forces the J — Ka 
colour to remain constant or to get bluer during light decline. 
The levitation of matter above the photosphere found here, 
along with the discovery of the development of a w arm chro¬ 
mosphere during rising light bv lWood et al.l ||2004|L provides 
strong evidence that the LSP phenomenon is associated with 
mass ejection from the photosphere of red giant stars. 

We examined two possible models for the explanation 
of the LSP phenomenon by comparing these models with 
broadband observations of both oxygen-rich and carbon 
stars. The first model assumes LSP light variations are due 
to dust absorption by ejected dust shells while the second 
model assumes a rotating star with a dark spot. 

In the dust shell ejection model for LSP stars, it was 
found that dust absorption by the ejected mass shell could 
not explain the observed variations in the (/,J — Ka) or 
{I,V — I) diagrams for either oxygen-rich or carbon stars. In 
oxygen-rich stars, the models fail to reproduce the variations 
in the (/,J — Ka) because the models only include variable 
dust absorption and not the variable H 2 O absorption noted 
above. The failure of the dust shell models to reproduce 
the observed variations in the {I,V — I) diagram for either 
oxygen-rich or carbon stars, and the failure to reproduce the 
observed variations in the {I,J — Ka) diagram for carbon 
stars, suggests that dust absorption in ejected mass shells is 
not the cause of LSP light variations. 

The LSP model involving a rotating, spotted star was 
not able to reproduce the observed variations in oxygen-rich 
stars in either the (/,/ — J) or {I,V — I) diagrams (or the 
7, J — Ka diagram, but here the J — Ka colour is affected 
by H 2 O absorption that is not included in the models). We 
thus conclude that the rotating, spotted star model is not 
the explanation for LSP light variations either . This result 
strengthens the conclusion of lOlivier fc WoodI ll2003l') that 
the rotation periods of LSP stars are too long for a spotted 
star model to explain LSP variations. A simple combined 
dust and spotted star model was no better at explaining 
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Figure 14. The variations of I with V — I, I — J and J — K for the models listed in Table [5] Models for oxygen-rich stars are shown 
in the upper 2 panels and models for carbon-rich stars are shown in the lower 2 panels. Solid lines with filled squares represent the dust 
shell model alone while solid lines with filled triangles represent the dark spot model alone. The combined dust and rotating spotted 
star models are shown as coloured lines with the colour indicating the phase lags 0. Note that the coloured lines have been arbitrarily 
shifted in I to prevent confusing overlap since it is only the slope of the lines that is of interest. The slopes of the three solid black lines 
in each diagram represent the observed variations. The central line corresponds to the mean value of the observed line slope for stars 
with a relative error in the slope of less than 0.2 and the outer lines vary from the mean slope by 1 a. 


the observed magnitude-colour variations than the dust and 
spotted star models individually. 


Our conclusion is that some process other than dust 
shell ejection or spot rotation lies behind the light variations 
seen in stars with LSPs. This unknown process seems to be 
responsible for the ejection of matter near the beginning of 
light decline. The ejected matter may be the origin of the 
excess mid- infrared emission due to dust that is seen in stars 
with LSPs dWood fc Nichollj|2009l L 
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